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Contrary to earlier expectations, several cosmic ray events with energies above 10 20 eV have been reported by 
a number of ultra-high energy cosmic ray observatories. According to the AGASA experiment, the flux of such 
events is well above the predicted Greisen-Zatsepin-Kuzmin cutoff due to the pion production of extragalactic 
cosmic ray protons off the cosmic microwave background. In addition to the relatively high flux of events, the 
isotropic distribution of arrival directions and an indication of small scale clustering strongly challenge all models 
proposed to resolve this puzzle. We discuss how the GZK cutoff is modified by the local distribution of galaxies 
and how astrophysical proton sources with soft injection spectra are ruled out by AGASA data. Sources with hard 
injection spectrum are barely allowed by the observed spectrum. If the most recent claims by AGASA that the 
highest energy events are due to clustered nuclei are confirmed, the most plausible explanation are astrophysical 
sources with very hard spectra such as extragalactic unipolar inductors. In addition, extragalactic magnetic fields 
need to be well below the current nano-Gauss upper limits. Alternatively, if the primaries are not nuclei, the need 
for new physics explanations is paramount. We present an overview of the theoretical proposals along with their 
most general signatures to be tested by upcoming experiments. 



1. Introduction 

The future of ultra high energy cosmic ray 
physics looks extremely promising. The present 
state of observations is particularly puzzling and 
the necessary experiments to resolve these puz- 
zles will be operating in the very near future. 
The puzzles begin with the lack of the pre- 
dicted Greisen-Zatsepin-Kuzmin (GZK) cutoff [jl], 
Contrary to earlier expectations, cosmic rays 
with energies above 10 20 eV have been detected 
by a number of experiments (for a review see || 
and for a more recent update see Q). 

If these particles are protons, they are likely to 
originate in extragalactic sources, since at these 
high energies the Galactic magnetic field can- 
not confine protons in the Galaxy. However, 
extragalactic protons with energies above a few 
times 10 19 eV can produce pions through inter- 
actions with the cosmic microwave background 
(CMB) and consequently lose significant amounts 
of energy as they traverse intergalactic distances. 
Thus, in addition to the extraordinary energy re- 



quirements for astrophysical sources to acceler- 
ate protons to ^ 10 20 eV, the photopion thresh- 
old reaction suppresses the observable flux above 
~ 10 20 eV. These conditions were expected to 
cause a natural high-energy limit to the cosmic 
ray spectrum known as the GZK cutoff 

As shown by the most recent compilation of 
the AGASA data the spectrum of cosmic 
rays does not end at the expected GZK cut- 
off. The significant flux observed above 10 20 eV 
together with a nearly isotropic distribution of 
event arrival directions Q challenges astrophys- 
ically based explanations as well as new physics 
alternatives (see |7|^] and references therein). In 
addition, the reported small scale clustering ||] 
tends to rule out most scenarios. 

This challenging state of affairs is stimulating 
both for theoretical investigations as well as ex- 
perimental efforts. The explanation may hide in 
the experimental arena such as an over estimate 
of the flux at the highest energies. This expla- 
nation has been proposed by the HiRes collab- 
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oration based on an analysis of their monocular 
data H . Even if this were the case, events above 
the GZK cutoff are also observed by HiRes. The 
Mono HiRes data looks more like the GZK fea- 
ture, as discussed in the next section, followed 
by indications of new sources at energies above 
the feature. Events past 10 20 eV pose theoretical 
challenges which will be explained in the future by 
either astrophysically novel sources or new funda- 
mental physics. 

2. The GZK Feature 

Ultra-high energy cosmic rays (UHECRs) prop- 
agating from extragalactic sources to Earth suffer 
a number of losses in their interaction with cos- 
mic backgrounds. These include pair production, 
photopion production, and adiabatic losses due 
to the universe's expansion. The adiabatic and 
pair production losses change the spectrum con- 
tinuously, while the pion production process has 
a threshold given by the pion mass, i.e., protons 
need to be at ultra-high energies to be able to 
produce a pion off CMB photons. A careful cal- 
culation would show that the GZK cutoff is not 
an absolute end to the cosmic ray spectrum but 
it generates a clear feature around 5 x 10 19 eV set 
by the pion mass threshold. The spectrum recov- 
ers at energies above this feature and the local 
distribution of sources can significantly affect the 
agreement between predicted and observed spec- 
tra (see JT(]{Ll|). 

To see this effect more clearly, we developed 
a numerical code that can calculate the UHECR 
spectrum for a given input spectrum (assumed to 
be a power law at the source with spectral in- 
dex 7, J(E) (x E a given spatial distribution 
of sources, and a choice for redshift evolution of 
sources JyJ. Some results from this study are 
shown in Figs. 1,2,4 and 8. 

In principle, a local overdensity of sources can 
decrease the gap between observed and detected 
events above the GZK cutoff. This effect can eas- 
ily be understood: photopion energy losses limit 
the maximum distance at which sources can con- 
tribute at the highest energies to a few tens of 
Mpc, while cosmic rays below the pion produc- 
tion threshold come from much larger volumes. A 
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Figure 1. Simulated fluxes for the AG AS A statis- 
tics of 728 events above 10 19 eV, and 7 = 3, us- 
ing a homogeneous source distribution (\ hatches) 
and the PSCz distribution (dense / hatches). The 
solid and dashed lines are the results of the ana- 
lytical calculations for the same two cases. The 
dash-dotted and dash-dot-dot-dotted lines trace 
the mean simulated fluxes for the homogeneous 
and the PSCz cases. 



local overdensity will increase the observed flux at 
the highest energies relative to the lower energy 
flux. If the UHECR source distribution is propor- 
tional to that given by the galaxy distribution, 
the observed local overdensity is not high enough 
to explain the data 0. We reached that con- 
clusion using the source distribution as given by 
the CfA and PSCz galaxy redshift surveys (PSCz 
results are shown in Figs. 1 and 2) . The local den- 
sity is only about a factor of two above the mean 
which is not enough to bridge the gap, although 
a priori UHECR sources may cluster differently 
from luminous matter (see, e.g., []l2"|| ). 

In order to explain the highest energy events 
with extragalactic protons, a more profitable av- 
enue is to study the effect of the spectral index. In 
Fig. 1, we set the input spectral index 7 = 3 such 
that the events below 10 19 5 eV also fit the pre- 
dicted flux. If the highest energy events are con- 
firmed by larger experiments, the evidence points 
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Figure 2. Simulated fluxes for the AG AS A statis- 
tics of 728 events above 10 19 eV, and 7 = 2.f, 
using a homogeneous source distribution with 
Zmax = 0.1 (/ hatches), the PSCz distribution 
with z max = 0.1 (horizontal hatches), and a ho- 
mogeneous source distribution with z max = 1 (\ 
hatches) . 



towards a new source with a spectrum harder 
than J(E) oc E~ 3 . Given earlier AGASA data 
as of ref. |L3| with 728 events above 10 19 eV, 
Fig. 2 shows how a spectral index 7 = 2.1 can 
more easily fit the data at the highest energies. 
Fig. 2 shows both the mean and the fluctuations 
around the mean flux for the data set in |l3| . As 
more events are observed, such as in the Auger 
Project jL4|, the fluctuations about the mean de- 
crease as shown in Fig. 8. In Fig. 3, we show the 
mean behavior of the GZK feature for different 
input spectra assuming a homogeneous distribu- 
tion of sources. As the spectral index decreases, 
the agreement with the observations improves. In 
sum, sources of UHECRs distributed as ordinary 
galaxies are marginally consistent with present 
spectral data and, for hard injection spectra, the 
GZK cutoff is not really a cutoff but a feature in 
the high-energy cosmic ray spectrum. 

The transition from the lower energy (~ 10 19 
eV) region of the UHECR spectrum to the higher 
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Figure 3. Propagated spectrum for source spec- 
tral index of 7 = 1.5, 2.1, 2.7. 



energy region (~ 10 20 eV) can be made smoother 
if the sources have a redshift evolution. In Fig. 
4, we show the effect of a source luminosity evo- 
lution for sources with luminosity oc (1 + z) m . If 
the source luminosity increases with redshift z, 
(m > 0), the flux of UHECRs at energies below 
the GZK cutoff will increase relative to the flux 
above the cutoff. This effect changes the shape of 
the spectrum, broadening the transition region. 

In addition to the presence of events past the 
GZK cutoff, there has been no clear counterparts 
identified in the arrival direction of the highest 
energy events. If these events are protons or pho- 
tons, these observations should be astronomical, 
i.e., their arrival directions should be the angular 
position of sources. At these high energies the 
Galactic and extragalactic magnetic fields should 
not affect proton orbits significantly so that even 
protons would point back to their sources within a 
few degrees. Protons at 10 20 eV propagate mainly 
in straight lines as they traverse the Galaxy since 
their gyroradii are ~ 100 kpc in fiG fields which is 
typical in the Galactic disk. Extragalactic fields 
are expected to be <C /iG, and induce at most 
~ 1° deviation from the source Jig] . Even if 
the Local Supercluster has relatively strong fields, 
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Figure 4. Simulated fluxes for the AGASA statis- 
tics of 728 events above 10 19 eV, and 7 = 2.1, 
using a homogeneous source distribution with 
Zmax — 1 and to = (horizontal hatches), m = 2 
(/ hatches), and m = 4 (\ hatches). 



the highest energy events are expected to deviate 
at most ~ 10° jl(|[l7). At present, no correla- 
tions between arrival directions and plausible op- 
tical counterparts such as sources in the Galactic 
plane, the Local Group, or the Local Superclus- 
ter have been clearly identified. Ultra high energy 
cosmic ray data are consistent with an isotropic 
distribution of sources in sharp contrast to the 
anisotropic distribution of light within 50 Mpc 
from Earth. 

In addition to the overall isotropic distribu- 
tion of arrival directions, the AGASA data shows 
evidence of a small scale clustering of events. 
The number of observed double and triple events 
seems much higher than expected for a random 
distribution ||. These clusters do not correlate 
with any known nearby galaxy population and 
give most models the hardest hurdle to overcome. 
The only positive cross correlation found thus far 
is between the UHECR clusters and very distant 
BL Lacertae objects J[8j| . These highly active 
galaxies are likely to be prime accelerators, but 
their location is too far from Earth to avoid a 



GZK cutoff to the spectrum. This positive cross 
correlation has inspired new physics proposals 
|T^| that make use of new neutral particles |2(| 
that need to be very long lived to be able to tra- 
verse cosmological distances. A more radical op- 
tion would be the breaking of Lorentz invariance 
that may render the neutron stable |pl| . Both 
possibilities give future experiments the perfect 
carrot. 

3. Astrophysical Zevatrons 

The puzzle presented by the observations of 
cosmic rays above 10 20 eV have generated a num- 
ber of proposals that can be divided into As- 
trophysical Zevatrons and New Physics models. 
Astrophysical Zevatrons are also referred to as 
bottom-up models and involve searching for accel- 
eration sites in known astrophysical objects that 
can reach ZeV energies. New Physics proposals 
can be either hybrid or pure top-down models. 
First we discuss astrophysical Zevatrons in this 
section and new physics models in the next. 

Cosmic rays can be accelerated in astrophys- 
ical plasmas when large-scale macroscopic mo- 
tions, such as shocks, winds, and turbulent flows, 
are transferred to individual particles. The max- 
imum energy of accelerated particles, -Emax, can 
be estimated by requiring that the gyroradius of 
the particle be contained in the acceleration re- 
gion: E max — ZeBL, where Ze is the charge 
of the particle, B is the strength and L the co- 
herence length of the magnetic field embedded in 
the plasma. For -E max ^10 20 eV and Z ~ 1, the 
only known astrophysical sources with reasonable 
BL products are neutron stars, active galactic nu- 
clei (AGNs), radio lobes of AGNs, and clusters of 
galaxies. Fig. 5 (know as a Hillas plot ||^]) high- 
lights the B vs. L for these objects. 

Clusters of Galaxies: Cluster shocks, although 
very large, are not able to accelerate protons to 
energies above ~ 10 19 eV |23|| . Propagation in 
the cluster also generates a GZK feature. 

AGN Radio Lobes: Jets from the central black- 
hole of an active galaxy end at a termination 
shock where the interaction of the jet with the 
intergalactic medium forms radio lobes and 'hot 
spots'. Of special interest are the most power- 
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Figure 5. B vs. L, for E max = 10 20 eV, Z = 1 
(dashed line) and Z = 26 (solid line). 



ful AGNs where shocks can accelerate particles 
to energies well above ~ 10 18 eV via the first- 
order Fermi mechanism p4| , p5[ . A nearby spe- 
cially powerful source may be able to reach en- 
ergies past the cutoff and fit the observed spec- 
trumlQ. However, extremely powerful AGNs 
with radio lobes and hot spots are rare and far 
apart and are unlikely to match the observed ar- 
rival direction distribution. If M87 is the primary 
source of UHECRs a concentration of events in 
the direction of M87 should be seen. The next 
known nearby source after M87 is NGC315 which 
is already too far at a distance of ~ 80 Mpc. 
Any unknown source between M87 and NGC315 
would likely contribute a second hot spot, not the 
observed isotropic distribution. The very distant 
radio lobes will contribute a GZK cut spectrum 
which is also not observed. 

The possibility of stronger Galactic and extra- 
galactic magnetic fields may reduce the problem. 
In particular, a strong Galactic wind can sig- 
nificantly alter the paths of UHECRs such that 
the observed arrival directions of events above 
10 20 eV would trace back to the North Galac- 
tic Pole which is close to the Virgo cluster where 
M87 resides [p7| . The proposed wind would fo- 



cus most observed events within a very narrow 
energy range into the northern Galactic pole and 
render point source identification fruitless. Full 
sky coverage of future experiments will be a key 
discriminator of such proposals. 

AGN - Central Regions: The powerful engines 
that give rise to the observed jets and radio lobes 
are located in the central regions of active galax- 
ies and are powered by the accretion of matter 
onto supermassive black holes. The central en- 
gines might themselves be the UHECR acceler- 
ators pq] . The nuclei of generic active galaxies 
(not only the ones with radio lobes) can accel- 
erate particles via a unipolar inductor not unlike 
the one operating in pulsars. In the case of AGNs, 
the magnetic field may be provided by the in- 
falling matter and the spinning black hole horizon 
provides the imperfect conductor for the unipolar 
induction. This proposal has to face the debilitat- 
ing losses that UHE charged particles face in the 
acceleration region due to the intense radiation 
field present in AGNs. In addition, the spatial 
distribution of objects implies a GZK cutoff of 
the observed spectrum. This limitation due to 
energy losses has led to the proposal that quasar 
remnants, supermassive black holes in centers of 
inactive galaxies, are more effective UHECR ac- 
celerators |29|| . From Figure 1-3, these models can 
only succeed if the source spectrum is fairly hard 
(7^2). 

Neutron Stars: Neutron star not only have the 
ability to confine 10 20 eV protons, the rotation 
energy of young neutron stars is more than suffi- 
cient to match the observed UHECR fluxes | g0| . 
However, ambient magnetic and radiation fields 
induce significant losses inside a neutron star's 
light cylinder. However, the plasma that ex- 
pands beyond the light cylinder is free from the 
main loss processes and may be accelerated to 
ultra high energies. In particular, newly formed, 
rapidly rotating neutron stars may accelerate iron 
nuclei to UHEs through relativistic MHD winds 
beyond their light cylinders ^jj] . This mecha- 
nism naturally leads to very hard injection spec- 
tra (7 ~ 1). In this case, UHECRs originate 
mostly in the Galaxy and the arrival directions 
require that the primaries be heavier nuclei. De- 
pending on the structure of Galactic magnetic 
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fields, the trajectories of iron nuclei from Galac- 
tic neutron stars can be consistent with the ob- 
served arrival directions of the highest energy 
events (see, e.g., 32 33[|). This proposal should 
be constrained once the primary composition is 
clearly determined. 

Gamma-Ray Bursts: Transient high energy 
phenomena such as gamma-ray bursts (GRBs) 
may also be a source of ultra-high energies pro- 
tons |34|,|35| . In addition to both phenomena hav- 
ing unknown origins, GRBs and UHECRs have 
other similarities that may argue for a common 
source. Like UHECRs, GRBs are distributed 
isotropically in the sky, and the average rate of 
7-ray energy emitted by GRBs is comparable to 
the energy generation rate of UHECRs of energy 
> 10 19 eV in a redshift independent cosmologi- 
cal distribution of sources. However, recent GRB 
counterpart identifications argue for a strong cos- 
mological evolution for GRBs. The distribution 
of UHECR arrival directions and arrival times 
argues against the GRB-UHECR common ori- 
gin. Events past the GZK cutoff require that 
only GRBs from <50 Mpc contribute. Since less 
than about one burst is expected to have occurred 
within this region over a period of 100 yr, the 
unique source would appear as a concentration of 
UHECR events in a small part of the sky. In ad- 
dition, the signal would be very narrow in energy 
AE/E ~ 1. Again, a strong intergalactic mag- 
netic field can ease the arrival direction difficulty 
dispersing the events of a single burst but also 
decreasing the flux below the observed level. 

4. New Physics Models 

The UHECR puzzle has inspired a number of 
different models that involve physics beyond the 
standard model of particle physics. New Physics 
proposals can be top-down models or a hybrid of 
astrophysical Zevatrons with new particles. Top- 
down models involve the decay of very high mass 
relics that could have formed in the early uni- 
verse. 

The most economical among hybrid propos- 
als involves a familiar extension of the standard 
model, namely, neutrino masses. If some flavor 
of neutrinos have mass (~ 0.1 eV), the relic neu- 



trino background is a target for extremely high 
energy neutrinos to interact and generate other 
particles through the Z-pole |3^,^]. This pro- 
posal requires very luminous sources of extremely 
high energy neutrinos throughout the universe. 
Neutrino energies need to be <U0 21 eV which im- 
plies primary protons in the source with energies 
>10 22 eV. The decay products of the Z-pole in- 
teraction are dominated by photons, which gives 
a clear test to this proposal. In addition, the neu- 
trino background only clusters on large scales, so 
the arrival direction for events should be mostly 
isotropic. Preserving a small scale clustering may 
be another challenge to this proposal. 

If none of the astrophysical scenarios or the 
hybrid new physics models are able to explain 
present and future UHECR data, the alterna- 
tive is to consider top-down models. The idea 
behind these models is that relics of the very 
early universe, topological defects (TDs) or su- 
perheavy relic (SHR) particles, produced after or 
at the end of inflation, can decay today and gen- 
erate UHECRs. Defects, such as cosmic strings, 
domain walls, and magnetic monopoles, can be 
generated through the Kibble mechanism as sym- 
metries are broken with the expansion and cool- 
ing of the universe. Topologically stable defects 
can survive to the present and decompose into 
their constituent fields as they collapse, annihi- 
late, or reach critical current in the case of su- 
perconducting cosmic strings |J8|,[39). The decay 
products, superheavy gauge and higgs bosons, de- 
cay into jets of hadrons, mostly pions. Pions in 
the jets subsequently decay into 7-rays, electrons, 
and neutrinos. Only a few percent of the hadrons 
are expected to be nucleons. Typical features 
of these scenarios are a predominant release of 
7-rays and neutrinos and a QCD fragmentation 
spectrum which is considerably harder than the 
case of Zevatron shock acceleration. 

ZeV energies are not a challenge for top-down 
models since symmetry breaking scales at the 
end of inflation typically are 3> 10 21 eV. Fit- 
ting the observed flux of UHECRs is harder since 
the typical distances between TDs is the Hori- 
zon scale or several Gpc. The low flux hurts 
proposals based on ordinary and superconduct- 
ing cosmic strings which are distributed through- 
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Figure 6. roton and 7-ray fluxes from necklaces 
for m x = 10 14 GeV (dashed lines), 10 15 GeV 
(dotted lines), and 10 16 GeV (solid lines) normal- 
ized to the observed data. 7-high and 7-low cor- 
respond to two extreme cases of 7-ray absorption. 



Figure 7. SHRs or monopolia decay fluxes (for 
mx — 10 14 GeV): nucleons from the halo (pro- 
tons), 7-rays from the halo (gammas) and extra- 
galactic protons. Solid, dotted and dashed curves 
correspond to different model parameters. 



out space. Monopoles usually suffer the opposite 
problem, they would in general be too numer- 
ous. Inflation succeeds in diluting the number 
density of monopoles and makes them too rare for 
UHECR production. Once two symmetry break- 
ing scales are invoked, a combination of horizon 
scales gives room to reasonable fluxes. This is 
the case of cosmic necklaces [Q which are hy- 
brid defects where each monopole is connected to 
two strings resembling beads on a cosmic string 
necklace. The UHECR flux which is ultimately 
generated by the annihilation of monopoles with 
antimonopoles trapped in the string In these 
scenarios, protons dominate the flux in the lower 
energy side of the GZK cutoff while photons tend 
to dominate at higher energies depending on the 
radio background (see Fig. 6). If future data 
can settle the composition of UHECRs from 0.01 
to 1 ZeV, these models can be well constrained. 
In addition to fitting the UHECR flux, topologi- 
cal defect models are constrained by limits from 
EGRET on the flux of photons from 10 MeV to 
100 GeV. 



Another interesting possibility is the proposal 
that UHECRs are produced by the decay of un- 
stable superheavy relics that live much longer 
than the age of the universe |^2j. SHRs may 
be produced at the end of inflation by non- 
thermal effects such as a varying gravitational 
field, parametric resonances during preheating, 
instant preheating, or the decay of topological 
defects. These models need to invoke special 
symmetries to insure unusually long lifetimes for 
SHRs and that a sufficiently small percentage de- 
cays today producing UHECRs. As in the topo- 
logical defects case, the decay of these relics also 
generates jets of hadrons. These particles behave 
like cold dark matter and could constitute a fair 
fraction of the halo of our Galaxy. Therefore, 
their halo decay products would not be limited 
by the GZK cutoff allowing for a large flux at 
UHEs (see Fig. 7). Similar signatures can oc- 
cur if topological defects are microscopic, such as 
monopolonia and vortons, and decay in the Halo 
of our Galaxy. In both cases the composition of 
the primary would be a good discriminant since 
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the decay products are usually dominated by pho- 
tons. In the case of SHR decays, the arrival di- 
rection distribution should be close to isotropic 
but show an asymmetry due to the position of the 
Earth in the Galactic Halo and the clustering 
due to small scale dark matter inhomogeneities 
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Figure 8. Simulated fluxes for the Auger pro- 
jected statistics of 9075 events above 10 19 eV, 
and 7 = 3, using a homogeneous source distri- 
bution (\ hatches) and the PSCz distribution (/ 
hatches). The solid and dashed lines are the re- 
sults of the analytical calculations for the same 
two cases. The dash-dotted and dash-dot-dot- 
dotted lines trace the mean simulated fluxes for 
the homogeneous and the PSCz cases. 



5. Conclusion 

Next generation experiments such as the Pierre 
Auger Project which is now under construction, 
the proposed Telescope Array, and the EUSO 
project and the OWL satellites will significantly 
improve the data at the extremely-high end of 
the cosmic ray spectrum. With these observato- 
ries a clear determination of the spectrum and 
spatial distribution of UHECR sources is within 



reach. In addition, the observations of UHE neu- 
trinos in horizontal showers promises to open a 
new window into the workings of our Universe. 

The lack of a GZK cutoff should become clear 
with Auger and most extragalactic Zevatrons may 
be ruled out. The observed spectrum will distin- 
guish Zevatrons from new physics models by test- 
ing the hardness of the spectrum and the effect of 
propagation. Fig. 8 shows how clearly Auger will 
test the spectrum independent of their clustering 
properties. The cosmography of sources should 
also become clear and able to discriminate be- 
tween plausible populations for UHECR sources. 
The correlation of arrival directions for events 
with energies above 10 20 eV with some known 
structure such as the Galaxy, the Galactic halo, 
the Local Group or the Local Supercluster would 
be key in differentiating between different models. 
For instance, a correlation with the Galactic cen- 
ter and disk should become apparent at extremely 
high energies for the case of young neutron star 
winds, while a correlation with the large scale 
galaxy distribution should become clear for the 
case of quasar remnants. If SHRs are responsi- 
ble for UHECR production, the arrival directions 
should correlate with the dark matter distribution 
and show the halo asymmetry. For these signa- 
tures to be tested, full sky coverage is essential. 
Finally, an excellent discriminator would be an 
unambiguous composition determination of the 
primaries. In general, Galactic disk models in- 
voke iron nuclei to be consistent with the isotropic 
distribution, extragalactic Zevatrons tend to fa- 
vor proton primaries, while photon primaries are 
more common for early universe relics. 

The hybrid detector of the Auger Project 
should help settle the present disparity between 
HiRes and AGASA by cross calibrating the two 
techniques. It will also determine the composition 
by measuring the depth of shower maximum and 
the ground footprint of the same shower. AGASA 
seems to detect a hint of composition shifts at the 
highest energies. This would be quite a surpris- 
ing development. In sum, the future looks very 
promising. The solution to the UHECR mystery 
as well as the birth of UHE neutrino astronomy is 
coming with the next generation of experiments 
which are under construction such as Auger or in 
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the planning stages such as the Telescope Array, 
EUSO, and OWL. 
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